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Abstract 

(<— ^ I An exact analytic calculation is presented for the spectrum of relic gravitational waves in the 
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^Sj ' scenario of accelerating Universe JIa + ^m = 1- The spectrum formula contains explicitly the 
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parameters of acceleration, inflation, reheating, and the (tensor/scalar) ratio, so that it can be 
employed for a variety of cosmological models. We find that the spectrum depends on the behavior 



OO 

^SJ ' of the present accelerating expansion. The amplitude of gravitational waves for the model Oa = 

'sj" ■ 0.65 is about ~ 50% greater than that of the model Oa = 0.7, an effect accessible to the designed 

^«0 , sensitivities of LIGO and LISA. The spectrum sensitively depends on the inflationary models with 

^J I a{T) oc |r|^^^, and a larger f3 yields a fiatter spectrum, producing more power. The current LIGO 

S:^ ' results rule out the infiationary models of /? > —1.8. The LIGO with its design sensitivity and 

-■"^ . the LISA will also be able to test the model of /3 = —1.9. We also examine the constraints on the 

»^| spectral energy density of relic gravitational waves. Both the LIGO bound and the nucleosynthesis 

O 

;— I ' bound point out that the model /? = —1.8 is ruled out, but the model /3 = —2.0 is still alive. 

Cd ■ The exact analytic results also confirm the approximate spectrum and the numerical one in our 



previous work. 

PACS numbers: 98.80.-k, 98.80.Es, 04.30.-w, 04.62- 



* Electronic address: |yzh@ustc. edu.cn 



I. INTRODUCTION 



Recently much progress has been made in the Laser Interferometer Gravitational waves 
Observatory (LIGO) with the typical sensitivity 10~^^ to 10~^^ being reached in the fre- 

n n n n 

quency range 100 ~ lOOOHz [l[ [2j [3| [J|. The chance to detect directly the gravitational 
waves (GW) has thus increased. Therefore, it is necessary to examine the possible objects 
of detections, such as the relic GW, which has a spectrum distributed over a rather broad 
range of frequencies. The stochastic background of relic GW has long been studied 5| [g] 
[3]. The calculations of spectrum generated during the transitions from the inflationary era 
;o the radiation-dominated era, or, to the matter-dominated era, have been carried out lal 
3 Q H H Q Q. Mo.e .centl^tud.. ha. been .ade on the effec. of the detailed 
slow-roll inflationary on the relic GW 15| [iGj [l7| , and on the other post- inflationary phys- 
ical effects on the relic GW [l8|. A constraint on the the tensor-to-scalar ratio r has been 
derived, using the CMB-galaxy cross-correlation [l9j. The relic GW can influence CMB and 
cause magnetic type of CMB polarizations, which can serve as another distinct signal of 



the relic GW. This kind of effects have been studied in Refs. 120 
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2^. On both 



iheoretical and observational issues of the relic GW, a recent review is given by Grishchuk 



25|. 



The observations on the SN la 26| 27|] indicate that the Universe is currently under 
accelerating expansion, which may be driven by the cosmic dark energy (^a ~ 0.7) plus 
the dark matter (fi^ ~ 0.3) with ^a + fi„ = 1 [281 ^ [Soi. The evolution of relic GW 
after being generated during the inflationary stage depends on the subsequent expansion 
behaviors of the spacetime background. The current accelerating expansion of Universe will 
have an impact on the relic GW and its spectrum. The spectrum of relic GW has been 
studied in specific models for dark energy, such as the Chaplyngin gas model 3l| and the 
X-fluid model [32]. In previous study we have studied the effects on the relic GW caused 
by the acceleration of the Universe for fixed Q\ = 0.7 and Qm = 0.3, and have obtained 



an approximate [33|, and a numerical spectrum [S^] of the relic GW. It was shown that, in 
comparison with the decelerating models, both the shape and amplitude of the spectrum 
have been modified due to the current accelerating expansion. However, in the previous 
work, the dependence of the spectrum upon the dark energy fraction ^a has not been 
examined. Extending these previous studies, in this paper we present an exact analytic 



calculation of the spectrum for any fraction Q\ of the dark energy. We will demonstrate 
how i^A affects the spectrum, discuss the dependence of spectrum upon the inflationary 
models. We will also examine the resulting spectrum by comparing with the sensitivity 
curves of the gravitational wave detections, such as the LIGO and LISA, and constrain the 
corresponding spectral energy density by the resent LIGO bound and by the nucleosynthesis 
bound. The resulting formula of spectrum will contain explicitly the parameter for the dark 
energy, as well as the parameters for the inflationary expansion, the reheating, the initial 
normalization of the amplitude, and the ratio of (tensor/scalar), so that it can be quite 
general and can be used in other possible applications. In this way the paper is also to 
serve as a useful compilation. Thus we have also listed the main formulae and the relevant 
specifications involved in the calculation of the spectrum. Throughout the paper we adopt 



notations similar to that of 10| |33| for convenience 



II. EXPANSION STAGES OF THE UNIVERSE 

The overall expansion of the spatially flat Universe is described by the Robertson- Walker 
metric ds"^ = a'^{T)[dT'^ — 6ijdx'^dx^], where r is the conformal time. The scalar factor a{T) 
is given by the following for various stages. 

The initial stage (inflationary) 

a{T) = Io\t 1^+^^, -oo < r < n, (1) 

where 1 + /? < 0, and ri < 0. The special case of /3 = — 2 is the de Sitter expansion of 
inflation. 

The reheating stage 

a{T) = a,{T - TpY+f% n<T<T,. (2) 

This stage is introduced to allow a general reheating epoch [lO| [33| . 
The radiation-dominated stage 

a{T) =ae{T-Te), Ts < T < T2. (3) 

The matter-dominated stage 

a{T) = a^{T - TrnY, T2<T < Te, (4) 



where te is the time when the dark energy density pA is equal to the matter energy density 
Pm- The redshift ze at the time te is given hy 1 + ze = (tt^Y^'^- If the current values 
^A ~ 0.7 and Q™ ~ 0.3 are taken, then 1 + ze ^ 1.33. For Q\ ~ 0.65 and ^2^ ~ 0.25 , then 



l + ;zs~ 1.23 |33|. 



The accelerating stage (up to the present time th) 

a{T) =lH\T-Ta\''^, Te<T<Th, (5) 

where the parameter 7 = 1.0 is the de Sitter acceleration for Q\ = 1 and Qm = 0. For 
the realistic model with ^a = 0.7 and Qm = 0.3 at present, we have numerically solved the 
Friedman equation 

{^f = H\Q^ + Qma~') (6) 

where a' = da{T)/dT. The resulting a{T) is plotted in FigHJ We have found that the 
expression of ([5]) with 7 = 1.05 gives a good fitting to the numerical solution a^r). Similar 
calculations show that 7 = 1.06 fits the model of ^a = 0.65 (in Figl2]), 7 = 1.048 fits 
the model ^a = 0.75 ( in FigJSD, and 7 = 1.042 fits the model ^a = 0.80. Thus, for the 
spatially fiat Universe (^a + ^m = 1), as long as the dark energy dominates over the matter 
component (^a > ^m), the generic fitting formula Eq.([5]) is effectively valid, and the range 
of values for the parameter 7 are close to 1.0. The constant r^ in Eq.(l5]) can be taken to be 
the same value, not very sensitive to the various values of ^a and Qm- 

There are ten constants in the above expressions of a(r), except (3, f3s, and 7 , that are 
imposed upon as the model parameters. By the continuity conditions of a{T) and ^(t)' at 
the four given joining points ri, r^, T2, and te, one can fix only eight constants. The other 
two constants can be fixed by the overall normalization of a and by the observed Hubble 
constant as the expansion rate. Specifically, we put q^th) = Ih as the normalization, i.e. 

\th -Tal = 1, (7) 

and the constant Ih is fixed by the following calculation 

As we have shown that 7 ~ 1.0 in the realistic models of acceleration expansion, so Ih is 
just the Hubble radius at present. Then everything in the expressions of a{T) from Eq.(l) 



through Eq.(5) is fixed up. For instance, one obtains 

^0 = iHhCE '' C2"^Cf Cl^"'^ , (9) 

where h=\l + (3\~^'+P\ (e = Te/th, C2 = {rE/nf, C = rs/r,, and Ci = {rJnf+P'. 

To completely fix the joining conditions we need to specify the time instants ri, r2, r^, and 
te that separate two consecutive expansion stages. From the consideration of physics of the 
Universe, we take the following specifications [33]: a{TH) / a{TE) = 1.33, a{TE) / a{T2) = 3454, 
a{T2)/a{Ts) = 10^^, and a(rs)/a(ri) = 300. From these, one makes use of the continuity 
conditions of a and a', and obtains 

\rE-ra\ = {1 + ze)~, 
2(1 + ze) 
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Tp 
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The above expressions all depend on the model parameters (3, jSg, and 7 explicitly, thus 

depend on Q;^. So we can expect that the spectrum of relic GW will depend on the present 

acceleration behavior of the Universe through 7. 

In the expanding Robertson- Walker spacetime the physical wavelength A is related to the 

comoving wave number k by 

A^^. (11) 

By Eq.(I7j) the wave number corresponding to the present Hubble radius is kn = 
2TTa{TH)/lH = 2tt. There is another wave number, kE = 27ra(r^)-?7 = kii/{l + ze)-, whose 
corresponding wavelength is the Hubble radius 1/ H at the time te- 



III. EQUATION OF GRAVITATIONAL WAVES 

Incorporating the perturbations to the Robert son- Walker metric, one writes 

ds'^ = a\T)[dT'^ - {6,j + hij)dx'dx% (12) 

where hij is 3 x 3 symmetric, representing the perturbations. The gravitational wave field is 
the tensorial portion of /ij^, which is transverse-traceless diK'^ = 0, 6'^^hij = 0, and the wave 
equation is 

d^{^/^d>'h,,{^,T)) = 0. (13) 

For a fixed wave vector k and a fixed polarization state a, the wave equation reduces to the 



second-order ordinary differential equation 33| 35 1 



a' 



C + 2-/^r + kX^' = 0, (14) 

where the prime denotes d/dr. Since the equation of /ijf (r) for each polarization a is the 
same, we denote h^ (t) by /ik(T) in the following. Once the mode function hk{T) is known, 
the spectrum h{k, r) of relic GW is given by 

h{k,T) = ^k\h,iT)\, (15) 

'IT 



which is defined by the following equation 

1"°° dk 

I h\k,T)—=<0\h'^ix,T)h,,i^,T)\0>, (16) 

where the right-hand-side is the vacuum expectation value of the operator h^^hij. The 
spectral energy density parameter ^g{k) of the GW is defined through the relation 



^/"<' 



where pg = 32^^ii,o^*fo i^ ^^^ energy density of the GW, and pc is the critical energy 
density. Then, one reads 

ng{k) = ^h\k,TH){^)\ (17) 

which is dimensionless. Note that the there might be divergences in the integration for p^, 
either infrared or ultraviolet. As is known, the infrared divergence is avoided if a infrared 
cutoff is introduced. This can be done since the very long waves with wavelengths comparable 



to, or longer than, the Hubble length do not contribute to the GW energy density 36|] . As for 



the very short wavelength portion, the ultraviolet divergences is also avoided by considering 



the Parker's adiabatic theorem [37|], which states that, during a transition between expansion 
epochs with a characteristic time duration At, the gravitons created will be suppressed for 
wavenumbers k > 1/ At. Thus, the spectrum segments in both the very low and very high 
frequency ranges should be discarded from these physical considerations. 

IV. INITIAL AMPLITUDE OF SPECTRUM 

Regarding to the relic GW, the initial conditions are taken to be during the inflationary 
stage. For a given wave number fc, the corresponding wave crossed over the horizon at a time 
Tj, i.e. when the wave length was equal to the Hubble radius: Aj = 27ra(rj)/A; to l/if(rj). 
From Eq.([T]) yields Hiji) = Iq^\1 + P\ ■ \Ti\^~^'^, and, for the case of exact de Sitter expansion 
of /3 = —2, one has Hiji) = l^^. Thus a different k corresponds to a different time Tj. Now 
choose the initial condition of the mode function hk{T) as 

l^'=(^^)l = 777T- (18) 



Then the initial amplitude of the spectrum is [10| [33| 

h{Kn) = A{^f^^ (19) 

where the constant 

A = 8v^6^. (20) 

The power spectrum for the primordial perturbations of energy density is P{k) oc \h{k, Th)\'^, 
and its spectral index n is defined as P{k) oc /c"^^. Thus one reads off the relation n = 2(3+5. 
The exact de Sitter expansion oi (3 = —2 leads to n = 1, yielding an initial spectrum 
independent of k, called the scale-invariant primordial spectrum. Other values of f3 will 
differ from the scale-invariant one. 

As is known, any calculation of the spectrum of the relic GW always has some overall 
uncertainty, originating from the normalization of the amplitude. Currently, from the obser- 
vational perspective, the best that one can do is to use the CMB anisotropics to constrain 
the amplitude, as they receive the contributions from both the scalar (density) and the ten- 
sorial (GW) primordial perturbations. However, there is a well known problem of how much 



relative contribution is from the relic GW, in comparison with the scalar type contribution 
(the density perturbations). There have been a number of discussion on the ratio of the relic 
GW to the scalar contribution, 

r = Ph/Ps. (21) 

Theoretically, it is, in our view, a problem of initial conditions on the ratio of the scalar 
and tensorial modes of comic perturbations. So far, in regards to the very long wavelength, 
some preliminary conclusion on the upper limit of GW contributions has been given, based 
upon the analysis on WMAP and the observational results of SDSS, for instance, r < 0.37 



(95% c.l.) [38(1 J39j. The final conclusion on this issue might be eventually rely on the 
more observations of CMB anisotropics and polarization (such as the Planck project in near 
future). In the following, the ratio r is treated as a parameter, representing the relative 
contribution by the relic GW to the CMB anisotropics AT/T at low multipoles. This 



will determine the overall factor A in flT9|) . Using the observed CMB anisotropics 29| is 
AT/T ~ 0.37 X 10^^ at / ~ 2, which corresponds to the anisotropics on the scale of Hubble 
radius, we put 

h{kH, th) = 0.37 X lO^V (22) 

Then the spectrum h{k, th) at the present time th is fixed. If we take the upper limit 
r = 0.37, then h{kH,TH) — 0.14 x 10^^. For smaller r, our calculation is still similar except 
the resulting spectrum is reduced by the corresponding numerical factor. 

V. ANALYTIC SOLUTION 

Writing the mode function hk{T) = jik{T)/a{T) in Eq. (fT^ . the equation for jJikij) becomes 

^^l + (A;2 - ^)/i, = 0. (23) 

For a scale factor of power-law form a{T) oc r", the general exact solution is of the following 
form 

/ifc(r) = Ci(A;r)5 4_i (fcr) + C2(A;r)2 Ji_^(A;r), 

where the constant Ci and C2 are to determined by continuity of the function ^kiT) and the 
time derivative {^ki'T') / o-i'^))' at the time instance joining two consecutive stages. 



The inflationary stage has the solution 

/ifc(r) = x2[yliJ^+i(a;) + v42J„(^+i)(a;)], -oo < r < ri, (24) 

where x = kr, and the two constants Ai and A2, determining the initial states, are taken to 
be 

A = V^/f^^'^'^'' ^2=^Ae-^, (25) 

cos pvr V 2 

both are independent of k. With Eg. (1251) the mode function /ifc(r) is proportional to the 

Hankel's function H \, 

fikir) = A^e-'^'' sm{P7i + '^)x"^ Hf^,{x) , (26) 

which, in the high frequency limit, approaches to the positive frequency mode 

lim //fc(r) -^ e-'^\ 

fc— >oo 

Thus the initial state fixed by Eq. (l25j) corresponds to the so-called adiabatic vacuum in the 



high frequency limit 401] 41 1. 



The reheating stage has 

^^k{T) = t--[B,Jp^^iit) + B2J_p^_iit)], n<T< r„ (27) 

where the variable t = kij — Tp), and the two coefficients Bi and B2 are fixed by joining 
the functions fik{T) and {nkij) / aij))' continuously at the time ti when the reheating epoch 
begins: 

^ _ [^ Vi(a^i)^/3,+|(ti)- V|(xi)J^^+i(ti) ^^ 

V ^1 '^/3. + i(^l)-^-/3,-|(^l) + '^-/3,-i(^l)^/3, + |(^l) 

V *1 Jp, + \itl)J-p,-litl) + ^-/3,-i(il)^/3, + |(^l) 

with Xi = kri, ti = k{Ti — Tp), and (1 + (3s)xi = (1 + (3)ti, which follows from the continuity 
of a{T) and a'{r) at the time Ti. 



The radiation-dominated stage has 

/ife(r) = Cle-*^ + C2e*^ r, < r < T2, (30) 

where the variable y = k{T — r^), and Ci and C2 are given by 



1 






Cl = ^^ { [(^ - -)-Jps + M + ■Jps^its)]Bl + [(^ - -)J„^, l(t.) - J_^, 3(t,)]i?2 



1 



piy2 _|_ p~3iy2 A p-W2 I p3iy2 



se \ ' ' ^ " S£ 

se 



(31) 
(32) 



C2 = -^-^ [-(^ + -)J,^^{ts) + J,^^{ts)]B, + [-{t + -)J-,.-|(t.) - J_,.-|(ts)]52 1^ , 

^^ L ys ys 

where t^ = k{Ts - r^), y^ = A;(rs - Te), and t^ = (1 + /JJy^. 
The matter-dominated stage has 

^k{T) = J — [D^Jl^{z)+D2J_:^{z)l T2<T<TE, (33) 

where z = k{r — Tm), and Di and -D2 are given by 



D, = [-e^y^ - —e^y^ + — ^^ ]Ci + [-e-*^^ + -^e"*^^ + ^j—]C2: (34) 



pW2 j p-W2 j 

2^/2 8|/| 2y2 8yi 

with 1/2 = ^^^2 - Te). 

The accelerating stage has 

TTS 

f^k{r) = ^ — [E^J^^is) + i?2J„^_i (s)], TE<r< th, (36) 

where s = k{T — r^), and Ei and E2 are given by 

i?i = A-i^ I j3(^^)[-:^-:dM _ j_^_3(s^)] - ^{zj,)j_^,.{se)\d, 

Se \ Se ^ ^ J 

+ { J^-^{^e)[ '-f^ J^,-i{se)\ + J_|(^i.) J_,_i(s^) [ D2, (37) 



i?2 = A-i^ <i J.{ze)J,^^^{se) - J|h^(s^) + VKs^)] J> Di 



+ <; -J.i{ze)J,^{se) - J_|[-^(se) + V|(si.)] ) D2. (38) 
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A = J,HisE)[ '-^ J-,-i{se)] - J-,-i{se)[-^^ + J,^{SE)] (39) 

'' Se Se 

where ze = k{TE - r^), se = Hte - r^), and 'JZe = -2s e- 

With all these coefficients having been fixed, the mode function hk^Tu) is known as a 
function of the wave number k at present time th-, so is the spectrum 

h{k,TH) = ^k\hk{TH)l (40) 

\/T[ 



as defined in Eq.( TT5l) . The above results form a useful compilation for computing the relic 
GW. To make use of the formulation fHOl) . one substitutes hkiru) = fi'k{'TH)/(^{'TH), where 
l^kiju) is given in Eq. (!36|) . Of course, to specify Hkiju), all the coefficients Ei , E2 through- 
out Ai, A2 have to be employed. One may, in his own computation, choose proper values 
of the parameters /3, /3s, and 7 for the specific expansion behavior, as well as the initial 
amplitude A in Eq.(l22]). 

For illustrations, taking the (tensor/scalr) ratio in Eq. fj2Tl) r = 0.37, we have plotted the 
exact spectrum h{k,TH) as a function of the frequency u = k/27ca in Fig |4] for 7 = 1.05, 
and in Figl5]for 7 = 1.06. In each of these figures of fixed 7, three spectra are shown for 
three infiationary models with j3 = —1.8, —1.9, and —2.0, and the parameter j3s = 0.598, 
—0.552, and —0.689 are taken, respectively (33|. As these figures show, the spectrum is 
scale-invariant with a fiat segment in the range u < lO^^^Hz and a slope segment in the 
range u > lO^^^Hz. 



Now we make a comparison o: 



the exact spectrum /i(z/, th) with the sensitivity curve 



from the recent S2 of LIGO Pj [2| [^ with the sensitivity 10^^^ to 10^^^ in the frequency 
range z/ = 10^ ~ lO^Hz. h{i/,TH) is given in Fig|6]for 7 = 1.05 and in Figl7]for 7 = 1.06. 
Both figures have plotted three spectra for infiationary models (3 = —1.8, (3 = —1.9, and 
/3 = —2.0, respectively. It is found that the infiationary models with /3 > —1.8 has an 
amplitude about an order higher than the LIGO sensitive curve. Even if we take a much 
lower value for the (tensor/scalar) ratio, say r = 0.05, the spectrum is still within the region 
detectable by the LIGO. Thus, the infiationary model (3 = —1.8 generating the relic GW 
with r > 0.05 is ruled out by the LIGO null results. The models (3 < —1.9 are still alive by 
this test alone. Moreover, when LIGO reaches its design sensitivity ~ 10~^^ in the frequency 
range in the forthcoming runs, it will also be able to test the model oi (3 = —1.9. 

Figj8]for 7 = 1.05 and Fig|9]for 7 = 1.06 give a comparison of the exact spectrum /i(z/, Th) 



with the sensitivity curve from LISA the Next Generation 42] in the lower frequency range 

11 



V = 10~^ ~ lO^Hz. It is interesting to notice that, when the LISA, as being designed, runs 
in space in the near future, it will be able to examine directly not only the model f3 = —1.8 
but also the model (3 = —1.9. For the latter model, even if a much lower value of the ratio 
r = 0.05 is taken, the LISA will still be able to detect it. This will be an improvement over 
the LIGO detection on the earth. However, as the two figures show, the inflationary model 
P = —2.0 seems to be still difficult to detect by the LISA with the capability as presently 
designed. 

Let us examine the dependence of the spectrum hiy^Tn) upon the dark energy f^A through 
the acceleration model parameter 7. In FigHni for a fixed /? = —2.0 we plot two spectra 
/i(z/, Tu) for the acceleration models 7 = 1.05 and 7 = 1.06 in a broad range of frequencies. 
As is seen, the difference between these two acceleration models are small. To show the 
details in the enlarged pictures, in Fig{T2] and Fig{TT] we have plotted the spectra in a 
narrow range of frequencies. It can be read that the amplitude of the model 7 = 1.06 is 
about ~ 50% greater than that of the model 7 = 1.05. That is, in the accelerating Universe 
with VL\ = 0.65 the amplitude of relic GW is ~ 50% higher than the one with ^a = 0.7. 
Note that the spectrum amplitude /i(z/, th) itself is very small, so this amount of ~ 50% of 
difference is probably difficult to detect at present. However, in principle, it does provide a 
new way to tell the dark energy fraction Q\ in the Universe. With the LIGO approaching 
its designed sensitivity, hopefully this difference can be detected thereby. As the LISA is 
currently designed, it will also be able detect this effect. 

Let us examine the spectral energy density fig(z^) and their constraints. FigHS] and 
FigHHare the plots of the spectral energy density fig(z^) defined in Eq. llTTl) for 7 = 1.05 and 
7 = 1.06, respectively. These plots of the exact analytic results agree with the numerical 
one in 34|. If we use the result LIGO third science run [Sj of the energy density bound for 
the fiat spectrum with Qq < 8.4 x 10^^ in the 69 — 156 Hz band, then the model /3 = —1.8 
is ruled out, but the models (3 < —1.9 survive. However, this LIGO constraint on the GW 
ener gy density is not as stringent as the constraint by the so-called nucleosynthesis bound 



431] [4^, whose main idea is the following: In the early Universe at a temperature T ~ a 
few Mev the nucleosynthesis process is going on. The relic GW will contribute to the total 
energy density p that drives the Universe expansion, thus will increase the effective number 
of species of particles (/*. More relic GW energy will enhance the freeze-out temperature 
for the process pe ^->- nu, and will lead to more neutrons available for the production of 

12 



heliuin-4 {^He). In practice the effective number of neutrino species N^ is used in place of 
g^.. Analysis has led to the nucleosynthesis bound on the relic GW energy density at the 
present time [431]: 

I Qg{u) d{log u) < 0.56 X 10~^ (41) 

where the value p^ ^ 2.481 x lO^^pc and the conservative value N^ < 4 have been used. 
Note that this is bound on the total GW energy density integrated over all frequencies. 
The integrand function should also have a bound Qgiiy) < 0.56 x 10~^ in the interval of 
frequencies 5(logi^) ~ 1. By this constraint it is also seen from FigJTSl and Fig{T4]that the 
model /3 = —1.8 has an Qglu) too high, is therefore ruled out, the same conclusion that we 
arrived at from FiglH] and FigH The model with jS = —1.9 are barely alive, as its energy 
density ^^(z^) tends to be growing higher with high frequencies. The model (3 = —2.0 are 
still robust since its spectral energy density fig(z^) is a flat function much lower than the 
limit in Eq.(|iT]). 



VI. ANALYTIC APPROXIMATION 



We now want to give an approximation to the above exact solution h{k, th) to recover the 



approximate analytic one given in 33j. The following approximation for the Bessel functions 
will be used 

M^) ~ V — cos(x - ^ - 7), a; > 1, (42) 

V TTX Z 4 

Note that the coefficients Di, D2, -Bi, B2, Ci, C2, Ei, and E2 are all functions of k, and they 
need to be approximated according to the value of k. 

In the long-wave limit, xi = kri <^ 1 and ti = (1 + /5s)a;i/(l + /5) <^ 1, from Eqs.(j28ll 
and (nU\i one has 

D, ^ 2-^+^=(i±|-)^+Hf-^=Ai, D2 ^ tf^^'^'A,. (44) 

D2 is higher order of ti and can be neglected in the following. 

From Eqs. (l31l) and fl32|) . in the long- wave limit t^ <^ 1 and y^ ^ 1? one has 

Si ^ itf^Di oc fc^, B2 ^ -Bi. (45) 
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From Eqs. fl34l) and (135|) , in the long-wave limit k ^ 1/t2, one has 

3? 
Ci ^ - — fii ~ A;^-\ C2«Ci, (46) 

SO 6*2 can be neglected. In the shortwave limit k ^ 1/t2, one has 

Ci ~ — 2ii?i sin2;2, C2 ~ 2zi?i cos2;2. (47) 

From Eqs. fl371) and fl38|) . for fcr^; <^ 1, one has 

Ei^Ci, E2^C2, (48) 

which also holds approximately for krE ^ 1, with some extra oscillating factors. 

With all these coefficients being estimated, now we can evaluate the approximation of 
the spectrum in Eq. flTBl) at the present time th, which is written as 

h{k,rn)=A^k\^'^^^^^ 



2Trb a^Tfj) 
Substituting the expressions Eq. (!36!) for fikiTn) and Eq.Q for Iq into the above leads to 



h{k,T^)=A[j{CE -" Cycler-] — \l^\EiJ,^{sn) + E,J_^.{sn)\. (49) 



fcHV 2 

Using the results from Eq. fl42l) through Eq. fHHl) . we approximate this expression by the 
leading term of power-law of k in various ranges of k. By some straightforward calculations, 
using I {th — Ta)/{TE — T2) | = 1/(1 + -2^), wc obtain the following expressions for the analytic 
approximate spectrum 

h{k,TH) = A{^f+P, k<kE; (50) 

h{k,TH) ^ A{-^f- \ ^ kE<k<kH; (51) 

h{k,TH)^A{^f J kH<k<k2; (52) 

,(,.,„).^(i.).«|_i_. ,,<,<,,; (53) 

where the small parameter e = (1 + /5)(1 — 7)/7, also depending on the behavior of the 
acceleration expansion through 7. The model 7 = 1 gives e = 0, and the results of Eqs. (l50|) 



through (|5^ reduce to exactly our early result given in [33|. The influence of detailed 
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accelerating expansion on the h{k,TH) is mainly demonstrated through the factor 1/(1 + 
ze)^~^'^, causing a variation in the magnitude of h{k, th)- For the inflationary expansion with 
(3^-2, the model of 7 = 1.05 (^a = 0.7) gives 1/(1 + ze)^^' = 0.423, and the model 
7 = 1.06 (r^A = 0.65) gives 1/(1 + ze)^'^'^ = 0.533, yielding roughly the amplitude of the 
model 7 = 1.06 greater than that of the 7 = 1.05 by about ~ 30%. The more accurate 
computation from the exact solutions shows an average difference of ~ 50%, as plotted in 
Figs.([l2]) and ([11]). Note that the factor 1/(1 + zeY = 0.987 for the model 7 = 1.05, and 
1/(1 + ZeY = 0.989 for the model 7 = 1.06, differing by only 0.2%, too small to tell by the 
current experimental detections. Therefore, in regards to the amplitude of relic GW, one 
can simply put e = in the approximate spectrum given in Eqs. (l50p through (I54p . just as 
it was in the model 7 = 1, causing only a difference of 0.2% in the amplitude for a variety 
of models with various 7. 

We remark that each of these expressions from Eq. (l5T]) to (15^ holds up to a numerical 
factor A, which contains certain oscillating factors of the form cos(A;r/f), or cos(y2) and 
sin(ts). In comparison with the decelerating models [lOj], Eq.( l5T]) is a new segment of spec- 
trum in kE < k < kn, whose occurrence is due to the acceleration of current expansion of 
the Universe. Besides, the three segments of spectrum, i.e., Eqs .f l52|) . fl53|) . and fl54l) . all 
have the extra factor (1 + ze)'^'"^ = (fim/^A)^^^^^ that are missing in the corresponding 
three segments in the decelerating models. 

VII. CONCLUSION 

We have presented a detailed calculation of the exact analytic spectrum of relic GW 
in the present fiat Q/^ + Q^ = 1 Universe in accelerating expansion. The resulting exact 
spectrum explicitly depends on the detailed behavior of the present accelerating expansion, 
characterized by the parameter 7 in the scale factor a{T) oc \t\^'^. It also explicitly depends 
on the inflationary model (3, the reheating model jSg, and the (tensor/scalar) ratio r as well. 
Therefore, the result is general enough to describe the GW spectrum /i(z/, th) produced 
from in a variety of accelerating cosmological models. One can use the formula in other 
applications by choosing a set of parameters (3, (3s, 7, and r. Besides, the analysis of the 
exact result gives the following conclusions: 

The GW amplitude of the model 7 = 1.06 is about ~ 50% greater than that of the 
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model 7 = 1.05, i.e., in the accelerating Universe with Qa = 0.65 the amplitude of relic GW 
is ~ 50% higher than the one with ^a = 0.7. Although it is probably difficult to detect 
at present, the effect does provide a new way to tell the dark energy fraction Q\ in the 
Universe. Hopefully this difference can be detected when the LIGO approaches its designed 
sensitivity ~ 10~^^, and the LISA runs in future. 

The spectrum depends sensitively on the parameter f3 of the inflationary models. A larger 
value of (3 yields a flatter spectrum /i(z/, th) with more power on the higher frequencies. 
The sensitivity curve of current LIGO rules out the inflationary models with (3 > —1.8. 
The LIGO with its design sensitivity and the LISA in future will also be able to test the 
/3 = — 1.9 model directly. 

The relic GW is also constrained through its spectral energy density fig(z^) by the resent 
LIGO bound and the nucleosynthesis bound. While both bounds rule out the inflationary 
model /3 = —1.8, the nucleosynthesis bound puts the model (3 = —1.9 in danger. However, 
the model /3 = —2.0 (de Sitter) is robust, since its spectral energy density ^g{i^) is flat and 
is ~ 10^^*^, much smaller than the nucleosynthesis bound. 

Finally, the exact analytic spectrum reduces to the approximate analytic and the numer- 
ical ones given in our previous study for the case 7 = 1. 
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FIG. 1: For the accelerating expansion with ^2^ = 0.7 the scale factor a^r) can be fitted by Eq.dS} 
with the parameter 7 = 1.05. 
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FIG. 2: For the accelerating expansion with JIa = 0.65 the scale factor a{T) can be fitted by Eq.dS]) 
with 7 = 1.06. 
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FIG. 3: For the accelerating expansion with il^ = 0.75 the scale factor aij) can be fitted by Eq.dS} 
with 7 = 1.048. 
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FIG. 4: For a fixed acceleration parameter 7 = 1.05 the exact spectrum hiy^Tn) is plotted for 
three inflationary models of /3 = —1.8, /3 = —1.9, and /? = —2.0, respectively. 
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FIG. 5: For a fixed acceleration parameter 7 = 1.06 the exact spectrum h[u,TH) is plotted for 
three inflationary models of /3 = —1.8, j3 = —1.9, and j3 = —2.0, respectively. 
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FIG. 6: For a fixed acceleration parameter 7 = 1.05 the exact spectrum h{u,TH) is plotted within 
the range of z/ = 10^ — lO^Hz for three inflationary models oi (3 = —1.8, (3 = —1.9, and (3 = —2.0, 
to compare with the sensitivity curve of second run from LIGO LI 4|. 
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FIG. 7: For a fixed acceleration parameter 7 = 1.06 the exact spectrum /i(i/, r//) is plotted within 
the range of z/ = 10^ — lO'^Hz for three inflationary models of /? = —1.8, (3 = —1.9, and f3 = —2.0, 
to compare with the sensitivity curve of second run from LIGO LI 4|. 
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FIG. 8: For a fixed acceleration parameter 7 = 1.05 the exact spectrum /i(i/, r//) is plotted within 
the range of z^ = 10~^ — lO^Hz for three inflationary models of /3 = —1.8, f3 = —1.9, and /3 = —2.0, 
to compare with the sensitivity from LISA the Next Generation 42]. 
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FIG. 9: For a fixed acceleration parameter 7 = 1.06 the exact spectrum /i(i/, r//) is plotted within 
the range of z^ = 10~^ — lO^Hz for three inflationary models of /3 = —1.8, /3 = —1.9, and /3 = —2.0, 
to compare with the sensitivity from LISA the Next Generation 42]. 
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FIG. 10: For a fixed inflationary parameter /3 = —2.0 the spectrum h[u,TH) is plotted for different 
acceleration models of 7 = 1.05 and 7 = 1.06. The two spectra are quite close to each other, and 
the difference in amplitudes of hiy^TH) is quite small, and difficult to tell in this figure. 
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FIG. 11: This enlarged picture is a portion of FigJlOi in the range u = 1 — 10'^ Hz to show fine 
differences in the spectrum h{v^TH) for different acceleration models. Note that the amplitude of 
h{i',T}{) for the model 7 = 1.06 is about ~ 50% higher than that of model 7 = 1.05. But in the 
range u = 10^ — lO^Hz the amplitude is only about < 3 x 10^^^, not accessible to the current LIGO 
yet. 
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FIG. 12: For a fixed /3 = —1.9 this enlarged picture in the range z^ = 10 — lO^Hz shows fine 
differences in the spectrum /i(i/, r//) for different acceleration models. Again the amplitude of 
h{i',TH) for the model 7 = 1.06 is about ~ 50% higher than that of model 7 = 1.05. Now in the 
range v = 10^ — 3 x lO'^Hz the amplitude is about ~ 10~^^, accessible to the LIGO as it approaches 
its designed sensitivity 10~^^. 
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FIG. 13: For a fixed 7 = 1.05 the spectral energy density r2g(z^) is plotted for the models of 
P = —1.8, /? = —1.9, and /? = —2.0. Obviously, the inflationary model of /3 = —1.8 has an 
ilg(i/) increasing too fast with the frequency z^, thus is ruled out by the LIGO bound and the 
necleosynthesis bound. ^g{i^) in the model of /? = —1.9 is narrowly below the necleosynthesis 
bound, but since ^g{i^) increases also too fast with v so it will barely survive. The model of /? = —2.0 
has a flat spectral energy density with a value ~ lO"^'', much smaller than the necleosynthesis 
bound. Thus the model (3 = —2.0 is robust. 
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FIG. 14: This picture is similar to Figll3lbut for a fixed 7 = 1.06. The conclusions are also similar 
to Figim 
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